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Abstract 


We apply the Adomian decomposition method (ADM) to obtain a subop- 
timal control for linear time-varying systems with multiple state and control 
delays and with quadratic cost functional. In fact, the nonlinear two-point 
boundary value problem, derived from Pontryagin’s maximum principle, is 
solved by ADM. For the first time, we present here a convergence proof for 
ADM. In order to use the proposed method, a control design algorithm with 
low computational complexity is presented. Through the finite iterations of 
algorithm, a suboptimal control law is obtained for the linear time-varying 
multi-delay systems. Some illustrative examples are employed to demonstrate 
the accuracy and efficiency of the proposed methods. 


AMS(2010): 49N05; 93C05. 


Keywords: Multiple time-delay systems; Pontryagin’s maximum principle; 
Adomian decomposition method. 


1 Introduction 


Optimal control of time-delay systems is one of the most challenging math- 
ematical problems in control theory. Indeed, the presence of delay makes 
analysis and control design much more complicated. Delays frequently oc- 
cur in mechanics, physics, population dynamics, biological, chemical, elec- 
tronic, and transformation systems; see [14]. The theory and the application 
of optimal control for linear time-delay systems have been developed per- 
fectly. However, as for nonlinear systems, synthesis problems that are solved 
by classical control theory lead to difficult computations. It is well known 
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that the nonlinear optimal control time-delay systems can be reduced to a 
two-point boundary value problem (TPBVP) involving both delay and ad- 
vance delay terms, implementing Pontryagin’s maximum principle (PMP); 
see [15]. In general, this TPBVP cannot be solved exactly and most. re- 
searches have been devoted to finding an approximate solution, for non- 
linear TPBVP. We briefly review some recent papers that are relevant to 
the method developed in the current work for time-delay optimal control 
problem. An averaging approximations for time-delay optimal control prob- 
lems [5], the B-spline approximation scheme [8, 24], the PMP [9], variational 
iteration method (VIM) [25-28], a novel feedforward-feedback suboptimal 
control of linear time-delay systems [13], Haar wavelets approach [29], hybrid 
of block-pulse functions and orthonormal basis [7, 11, 19, 22,30], composite 
Chebyshev finite difference method [20], the Hamilton-Jacobi-Bellman equa- 
tion [6], a delay-dependent stability of neutral systems [4], an interior-point 
algorithm [34,35], and an embedding process that transfers the problem to a 
new optimal measure problem [17]. 

The topic of the Adomian decomposition method (ADM) has been rapidly 
growing in recent years. It was first proposed by Adomian [1,2]. In this 
method, the solution of functional equations is considered as the sum of an 
infinite series usually converging to the solution. A lot of research works 
have been conducted recently in applying this method to a class of linear and 
nonlinear partial differential equations; see [33]. The Adomian’s decompo- 
sition has many advantages: It does not require any kind of discretization, 
linearization, or perturbation of the variables and the equation, therefore it 
does not need any modification of the actual model that could change the 
solution; it is efficient on providing an approximate or even exact solution 
in a closed form, to linear and nonlinear problems and provides a fast and 
accurate convergent series, and therefore it is only necessary to calculate a 
few terms of the series in order to obtain a reliable approximate solution; the 
method depends only on the known function wo(t) and the algorithm is of 
simple implementation. The method, has been widely applied to solve non- 
linear problems, and different modifications are suggested to overcome the 
demerits arising in the solution procedure [3,32]. 

This paper concerns with a class of nonlinear quadratic optimal control 
problem with multi-delay systems. Applying the main ideas of the shooting 
method to the basic and also an ADM. By applying the necessary optimality 
conditions, we obtain iterative formulas for the ADM. By using the finite-step 
iteration of algorithm, we obtain a suboptimal control law. The convergence 
of the ADM is studied and for illustrating the effectiveness of these methods, 
some test problems are investigated. Four illustrative examples are given to 
demonstrate the simplicity and efficiency of the proposed method. 

The structure of this paper is arranged as follows: Section 2 is devoted to 
Pontryagin’s maximum principle used for solving linear time-varying multi- 
delay system. Section 3 is dedicated to the proposed design approach to 
solve a closed loop optimal control problem based on the ADM and conver- 
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gence of the method is demonstrated. Section 4 is devoted to the suboptimal 
control strategy and algorithm for the proposed method. In Section 5, the 
numerical examples are simulated to show the reasonableness of our theory 
and demonstrate the performance of our network. Finally, we end this paper 
with conclusions in Section 6. 


2 Problem statement and optimality conditions 


Consider the following linear time-varying multi-delay system: 


a(t) = A(t)a(t) + Ai(t)x(t — Tr) + B(t)u(t) + Bi(t)u(t — t.), 
a(t)=9(t), to—tTe <t< to, (1) 
u(t) = ¥(#), to -—Tu St < to, 


where x(t) € R” and u(t) € R™, are the state and control vectors, respec- 
tively; A(t), A,(t), B(t), and B,(t) are real, piecewise continuous matrices 
of appropriate dimensions defined on the appropriate intervals; ¢(t) and w(t) 
are specified initial functions; 7, and 7, are constant positive scalars. Here, it 
is assumed that system (1) is controllable and that 7, <7. Find the control 
signal u(t) that minimizes the cost functional: 


J = Lo" t)Qyelts) +3 | “(a7 QQWalt) + u"(HR()ult)) dt, (2) 


where, the matrix Qs € R”*” is symmetric positive semidefinite, Q(t) € 
R"*”" and R(t) € R™*™ are chosen to be positive semidefinite and positive 
definite matrices, respectively. 


The Hamiltonian function for the problem is 


Hu = 52 ()Q()a(t) + set R(A)u(t (3) 
+ AP (£)[A(t)x(t) + Ar(t)x(t — Te) + B(t)u(t) + Bi (t)u(t — Ty)], 


where A(t) € R” is the vector of the Lagrange multiplier. According to the 
necessary conditions for optimality, we can obtain the following nonlinear 
TPBVP [15,31]: 


A(t)e(t) + Ai (t)e(t — T2) — (Si (t) + So(t))A(Z) 
—S3(t)A(t+ Tu) — Sa(HA(E- Tu), to <t<tp—Tu, 
A(t)x(t) + Ai(t)x(t — Te) — S1(t)A(t) 

—Sy(t)A(t = Ti) ty — Ty < t < tf 


(4) 


and 
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—Q(t)a(t) — AT(t)A(t) 
AD =k GA Gee teh: Heth e, (5) 
—Q(t)x(t) — AT (E)A(t), tp —Te Kt <ty 


with initial conditions 


x(t) — p(t), to — Tx < t 
u(t)=V(t), to-—m<t 


where 
S1(é) = BOR @B, 
So(t) = By(t)R7'(t— ty) BY (A), 
S3(t) = BY)R*() BU (t+ Tu), 
Sa(t) = Bi(t)R- +t Tu) BT (t = Ty); 


x(t —T) is the time-delay term and A(t +7) is the time-advance term. Also, 
the optimal control law is obtained by 

—R-1(t)B™ (t)A(t) 
ur(t)= 9 -RU(E)BI(t+tu)At+ Tu), to<t<tp—Tu, (7) 
—R-()BT(t)A), tp —Tu <b < ty. 


The optimal can be implemented as a closed loop optimal if the co-state 
vector obtained consists of linear function of the states and a nonlinear term, 
which is the adjoint vector sequence, in the form 


A(t) = P(t)a(t) + g(t), Ate) = Qrx(ty), (8) 


where P(t) € R”*” is an unknown positive-semidefinite function matrix and 
g(t) € R” is the adjoint vector. 


Substituting (8) into equation (4) yields 


&(t) = [A(Q) — Sit) PM] 2(t) — Si)g(t) + Ar()a(t — te) + F(d), 
a(t)=¢(t), to-TStX<to, (9) 


where 


F(t) = 9 —Sa(t) [P(t — tu)a(t — Tu) + 9(t— Tu)], to St < ty — Tus 
—Sa4(t) [P(t — Tu)a(t— Tu) + g(t —Tu)], tp—Tu <t< ty. 
(10) 
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Computing the derivatives to the both sides with respect to t of equation (8), 
we have 


Mt) = P(t)a(t) + P(t)e(t)+ g(t), to<t<ty 
= [Pw + P(t)A(t) — P(t)Si(t)P(O)| x(t) — PO) Si (ge) 
+ P(t)Ai(t)a(t — T,) + P(t) F(t) + g(t). (11) 
Putting (8) into equation (5), we get 
—Q(é)x(t) — AT()P)x(t) — ATHg) 


Mt) = ¢ —AT (t+ Te)[P(t + Te)a(t + Te) +.9(t+t:)],to<t<ty—T2, (12) 
—Q(t)x(t) — AT (t)P(t)x(t) — AT (t)g(t), tp — Te <t < ty. 


Thus, from (11) and (12), we can obtain the following Riccati matrix 
differential equation: 


— P(t) = P()A(t) + ATH) Pt) — PH)Si(t) P(t) + Q(t), Plts) = Qs, (13) 
and the following adjoint vector differential equation: 


g(t) = —[A®) — HO)PO] 9(H-POAi(Ha(t—7)+G(), g(t¢) = 0, (14) 


AT (t+ Tx) [P(t + tx) x(t + Tx) + g(t + T2)] + P(t)S2(t) [P(Ha(t) + g(t)] 
—P(t)S3(t) [P(t + tu)a(t + Tu) + g(t + Tu)] 

G(t) = P(t)Sa(t) [P(t — Tu) x(t — Tu) + g(t — Tu)], to <t<ty—Te, 
P(t)S2(t) [P()a(t) + g(t)] + P(t)S3(t) [P(t + Tu)a(t + Tu) + g(t + Tu)] 
P(t)S4(t) [P(t — Tu) a(t — Tu) + g(t — Tu)], tp —Te <t<tyf—Tu, 
P(t)S4(t) [P(t — Tu)x(t — Tu) + g(t — Tu)], tp —Tu KEK tp. 


(15) 
Substituting (8) into (7) yields 
—Ro*(t) BT (t)[P@a(t) + 9] 
t 


u*(t) = 4 —RV(t) BYE (t+ tu)[P(t+ Tu) a(t + Tu) + g(t + Tu)], to St < ty — Tu, 
—K*()B* @|P@)a@)+ g@)|,t¢— mu St< tp. 


For the sake of simplicity, let us define the right hand sides of (9) and 
(14) as follows: 


filt,v,g) = [A(t) — Si(t)P(@)] a(t) — Si(t)g(t) + Ar (a(t — Te) + F(t), (16) 


falt,2,9) == —[AQ) — Si) PO]? g(t) — PArMa(t- 72) + G(), (17) 
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where F'(t) and G(t) are relations (10) and (15), respectively. 
Thus the TPBVP in (9) and (14) changes to 


= fo(t, x, 9), (18) 


Note that, relations (18) form a nonlinear TPBVP with time-varying coef- 
ficient involving both delay and advance terms. The exact solution of this 
problem is, in general, extremely difficult, if not impossible. In the next sec- 
tion, we propose another analytic approximate method based on ADM, for 
this purpose. 


3 Adomian decomposition method 


In order to illustrate the basic concepts of the ADM, we consider the following 
equation: 


L(u) + R(u) + M(u) = A(t), (19) 
where u(t) is an unknown function, £ is a linear operator, that is assumed 
to be invertible, R is another linear differential operator, (uw) represents 
the nonlinear terms, and h is a continuous function. Applying the inverse 
operator £~! to both sides of (19) and using the given conditions, we obtain 


u=f—L(R(u))-L7*(N(u)), (20) 


where the function f(t) represents the terms arising from integrating the 
function h(t) and using the initial condition. 
The standard Adomian method defines the solution u(t) of (19) as a series 


u(t) = S> un(t), (21) 
n=0 


where the components u,(t) are usually determined recurrently. Substituting 
this infinite series into (20) leads to 


Yo un(t) = FE) — £7 (Ro a) ae (wo, a) - (22) 
n=0 n=0 n=0 


The nonlinear term in (21) can be computed by substituting 


N(u) = $5 An(uo, ta, -- +; tn); (23) 
n=0 


Solving linear optimal control problems of the time-delayed systems by ... 171 
where A,, is the Adomian polynomials, which can be determined by 


1a 
Ay ia, Dove] ’ n=1,2,3,...,. (24) 
: G4 


n=0 


Now, substituting (23) into (22) leads to 


So un(t) = f(t) -— £74 @» a) =f? (> 4s] 
n=0 n=0 n=0 


Each term of series (21) is given by the recurrent relation 


Now, we briefly describe how to apply the ADM to systems (18). For this 
purpose, we use a shooting method like procedure combine with the ADM 
for solving TPBVP in (18). 


Remark. It is necessary to transform the boundary value problem into 
a initial value problem. Therefore, we must find a € R such that the condi- 


tion g(t¢) = 0 can be replaced by the condition g(t9) = a. Thus we rewrite 
the TPBVP (18) as follows: 


= fo(t, z, 9), (25) 


where a € R is an unknown parameter. This parameter will be identified 
after sufficient iterations of ADM. 


Based on the ADM, we seek the solution {2, g} as follows: 


N N 

z= lim y x = lim y 
N-oco mA g N-co Ins 

n=0 n=0 


and hence the recursive relationship is found as 


t= Ais, wD, 
GeiHL Aan, wen, 
x(to) = Xo, g(to) = a, 


with inverse £~1(.) = fi (at and 
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f(t, @n,9n) = Ye Am k=, 2, 


where Aj», are the Adomian polynomials and are calculated by 


An = a whe (45 Xe me Sod n| , n=0,1,2,.... (26) 


q=0 


Take the first n + 1 terms of the nth approximation of x and g as follows: 


i =fo+ Dial (Ais), (27) 


Un = 90+ dyer £(A2,i-1)- 


Find the sequences ®, = %) +-:-+ a, and WV, = go +--+ + gn such that 


(28) 


©, = 29 + L7"(filt, ®n—1, Un_1)), n= 1, 
Vv, = go t+ Lo" (fo(t, Oni, Un-1)), n= ale 


where xo(t) = x(to) = @o, go(t) = g(to) = 


Theorem 1. Assume that {37/9 ae(t)} and {Sop_o ge (t)} are the solution 
sequences produced by ADM formula (28), which converge, respectively, to 
x(t,a) and g(t,a), as n + oo. Then Z(t, a), g(t,a) are the exact solutions 
of (25). Accordingly, Z(t, @),g(t,@) are the exact solutions of (18) when @ 
is the real root of g(tr,a) = 0. 


Proof. Let consider TPBVP (25) as follows: 


(to) + Lo iG r,9)I, 
(to) + L~"|fo(t, 2, 9)I, 


where L= “( yand £-) = =f ( -)dt. We have 


t 
Ln(E) =| [Ai n—1(%0, %1,°°* >Um—1590;91,°°° »In-1)] ds, n= 1, (29) 


to 
t 
9n(t) =| [Aon—1(%0,%1,°°° >Un—15,90;91,°°° »9n-1)] ds, n= 1, (30) 
to 
ro(t) = x(to) = 20, go(t) = g(to) =a. 


Taking limits of both sides of (29) and (30) as n > co implies 
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lim tp(t) = lim Ay p—1(%0,%1,°°* , LR “++ Gk—1)| ds 

om x(t) m ee 1,k i( 0,41; »>Tk-1,90,91; > 9k 1) ? 


n 


n t 
im, 2 gute) — [ i, So dancer >Uk-1,90,915°°° | ds. 


Then 


a(t,a) = | [f(s,2(s,4),d(s,4))las 


to 


g(t, a) = / Lfo(s, £(s, a), G(s, a))|ds. 


to 


Differentiating both sides with respect to t yields 


Moreover, if t = to, then from (29) and (30), tn(to) = 0, gn(to) = 0 for every 
n> 1. Thus 


n n 


Y= ex(to) = 2o(to) = 20, — ge(to) = go(to) = a, 


k=0 k=0 


or equivalently, Z(to,a@) = 20, g(to,a@) = a. Hence, £(t,a@) and g(t,a) are 
the exact solutions of (25). In addition, they are the exact solutions of (18), 
only if the condition g(t) = 0 is satisfied. So, it is straightforward to choose 
the unknown parameter a € R” such that g(tr,a) = 0. Denoting this real 
root of g(tr,«) = 0 by @ completes the proof. 


4 Suboptimal control design strategy 


Consider the linear time-varying multi-delay system (1) with cost functional 
(2). Then, the Nth order suboptimal trajectory-control pair is obtained as 
follows: 


tes = Dino te(O), (31) 


w= >i 9G) 


and 
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= | (t)an(t) + gn (E)] 
—R-1@) BP (t+ m)[P(t + tu)en(t + Tu) (32) 
ee ai) | to <t<tf —Tu, 
—R-'()BT(t)[P@)en() + 9n(t)], tp—Tu St < ty. 


The integer N in (31) and (32) is generally determined according to a concrete 
control precision. Then, the following cost functional can be calculated: 


t 
In = FehltQren(ty) +5 f' (CROQMWen() + uh OROu(®) a (33) 


The Nth order in (31) and (32) has the desirable accuracy, if for given 
positive constant € > 0, the following condition holds jointly: 


aH Ze, (34) 


Jn 


If the tolerance error bound is chosen small enough, the Nth order suboptimal 
control law will be very close to u*(t), and thus, the value of cost functional 
in (33) and its optimal value J* will be almost identical. 


Algorithm: Suboptimal control law of system (1): 


Step 1: Obtain P(t) from (13). Let xo(t) = x(to) = d(t), go(t) = g(to) 
and k = 1. 


Step 2: Compute x;(t) and g,(t) from (29) and (30). 
Step 3: Let N =k and obtain y(t) and un(t) from (31) and (382). 
oe In — In-1 
Step 4: Calculate Jy according to (33). If a <¢€, then stop 
N 


and output uy(t), go to step 5; else, replace k by & + 1 and go to step 2. 


Step 5: Stop the algorithm; xy(t) and uy(t) are accurate enough. 


5 Numerical examples 


In this section, the proposed method is illustrated by some test problems. 
The calculations are performed by using MATLAB software. 


Example 1. Consider the following linear time-varying multi-delay systems 
[20]: 
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_ 1 3 
&(t) =2(t— 5) +ta(t— 3) + u(Z), 0<t<l, (35) 
z(t) =t+1, —-2<t<0, 
with the cost functional 
a Ls 
J=-a°(1)4+= u* (t)dt (36) 
2 2 Jo 
The exact solutions of x(t) and u(t) are given by 
239075 43 , 3129081,2 1178769 1 
720332 ¢ / 1681328 ¢ 611392 t+, 0<t< 4? 
14,3 , 1119201,2 680513, , 7039811 1 1 
3° + -gr06e4 & 420332 + 7336704" 4St<3, 
— J 239075 44 | 337595943 2093255542 , 1156163, , 56216927 1 3 
x(t) iesi328¢ + 5oa3084¢ 13450624¢ + qo00784¢ + tera074s87 3 Sb <7, 
47815 45 4 2306773 44 _ 246187143 , 1486537742 4 17419475 4 
420332” — 10087968 2521992 53802496" ' 26901248 
2652913 3 
+774673408" a<t<l, 
and 
29689342 , 2078251 1484465 1 
720332 + 840664 t 611392 ? 0<t< 4? 
_ J 296893 890679 1 i 
u(t) 10166 — 4203327 a<t<9, 
296893 1 
~~ 210166? 3<t<l 


Furthermore, the exact value of the cost functional J is given by J = 
1.70648554. 

In order to obtain an accurate enough suboptimal control law, we applied the 
proposed algorithm with the tolerance bounds e = 2x 107°. In this case, con- 


io = J9} = 1757 19-6 < 
10 


vergence is achieved after 10 iterations, that is, 


2x 10-*. 


To analysis the accuracy and effectiveness of the ADM and VIM, the rel- 
ative error of objective value is summarized in Table 1, for several iterations. 
Accordingly, results show that in N = 10, ADM converges to the exact solu- 
tion. The optimal value of the cost functional is obtained as Jig = 1.706485. 


We compared the results obtained from ADM and VIM with exact solu- 
tion. Figure 1 shows the simulation curves of u(t) at iteration time 10 and 
the corresponding state trajectory x(t). Results of both methods are very 
close to exact solution as shown by Figure 1. This confirms that the proposed 
method yields excellent results. 


Example 2. Consider the following linear time-varying multi-delay systems 
[7, 16, 21, 36]: 
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Table 1: Simulation results of Example 1 at different iteration times. 
method ADM ADM VIM VIM 
Jn — JN_-1 Jn — Jn-1 
N J. ———_— J, ——____—_— 
N = N Te 


6 1.682312 — 1.678321 — 

7 1.691205 | 5.255 x 10~? | 1.685342 | 4.165 x 10~? 
8 1.703156 | 7.016 x 10~? | 1.692335 | 4.132 x 1073 
9 1.706484 | 1.961 x 10~? | 1.697346 | 2.952 x 1073 
10 1.706485 | 1.757 x 10~® | 1.697349 | 1.767 x 10~° 


K ADM 
oO VIM 
Exact 


1 1 2 
£(t) = —2(t) + a(t — 3) + ut) — gut — 3) 0<t<l, 
x(t) = 1, 5 <t<0, (37) 
uf) =0, —Z<t<0, 
with the cost functional 
Ie ari 15 
J== x(t) + =u“ (t) } (t)dt. (38) 
2 I, 2 


According to system (1), we have A = —1,A; = B=1,B, = -3,Q — 
i,-R= 5. 

In order to obtain an accurate enough suboptimal control law, we applied 
the proposed algorithm with the tolerance bounds € = 8.1 x 107". Simulation 
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results at different iteration times are summarized in Table 2. From Table 
2, it is observed that convergence is achieved after 10 iterations, that is, 
a" 49) — g. 040 x 10-7 < 8.1 x 1077. 

Jio 

In Table 3, the minimums of J using the hybrid of block pulse and Leg- 
endre polynomials [21], hybrid of general block pulse and Legendre polyno- 
mials [36], orthogonal basis [16], hybrid of block pulse and orthogonal Taylor 
series [7], and present two methods are listed. Also the suboptimal control 
and state of the proposed ADM and VIM are demonstrated in Figure 2. 


Table 2: Simulation results of Example 2 at different iteration times 


method ADM ADM VIM VIM 
Jn — JN-1 Jn — In-1 
N J —_————_ J — 
N Te N Te 


6 0.36430512 = 0.35935108 = 

7 0.36512980 | 2.258 x 1072 | 0.36431299 | 1.382 x 107? 
8 0.37311212 | 2.139 x 107? | 0.37511283 | 2.879 x 1072 
9 0.37311294 | 2.519 x 10-® | 0.37311210 | 5.362 x 1073 
10 0.37311291 | 8.040 x 10-7 | 0.37311310 | 2.680 x 10-6 


Table 3: The cost functional values for Example 2 


Method Cost functional values 

Marzban and Razzaghi [21] 0.37311241 
Wang [36] 0.37312682 

Kellat [16] 0.3731123 
Dadkhah and Farahi [7] 0.373112935 

Proposed method (N = 10) 

ADM 0.37311291 

VIM 0.37311310 


Example 3. Consider the following linear time-varying multi-delay systems 
[12, 18, 23]: 


£1 (t) = to(t) + 21 (t— 1), 7 
£(t) = tai(t ) + 2a1(t — 

x1(t) = r9(t) = i —= 
u(t)=5(t+1), 0. 


oe BH 


with the cost functional 
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Figure 2: The suboptimal control and state, when N = 10 for Example 2. 


u(t) 
(t + 2) 


J= 578) + 03(3) +2 | *202(2) + 2a (t)x2(t) + 23(t) + Jdt. (40) 
0 


2 


According to system (1), we have 


and from (2), we get 


a;=(q9): a=(3;). R=1/(t+2). 


In order to obtain an accurate enough suboptimal control law, we applied 
the proposed algorithm with the tolerance bounds e« = 2.3 x 1076. Simu- 
lation results at different iteration times are summarized in Table 4. From 
Table 4, it is observed that convergence is achieved after 12 iterations, that 
eee a 7a Jit) _ 9.970 x 10-6 < 2.3 x 10-8. Table 5, a comparison is made 

12 
between the value of J obtained by the present method with N = 12, to- 
gether with the value of J reported in the literature: Malek-Zavarei [18], by 
employing an iterative approach for determining the suboptimal control of 
the mentioned problem; the method of Hwang and Chen [12], by employing 
Legendre polynomials of order 20; and the method of Marzban and Pirmora- 
dian [23], by employing a direct approach based on a hybrid of block-pulse 
functions and Lagrange interpolating polynomials. 


is, 
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Table 4: Simulation results of Example 3 at different iteration times 


method ADM ADM VIM VIM 
Jn —JN_1 Jn — Jn-1 
N J J oe 
N i N Tig 
8 23.21056 - 22.17650 - 


9 23.05198 | 6.879 x 10-7 | 22.15431 | 1.001 x 10~? 
10 22.02108 | 1.681 x 10-3 | 22.01430 | 5.903 x 10-3 
11 22.02230 | 5.539 x 10-° | 22.02143 | 3.237 x 10-4 
12 22.02235 | 2.270% 10-° | 22:02115.| 1.271 x 10-* 


Table 5: The cost functional values for Example 3 


Method Cost functional values 
Malek-Zavarei [18] 24.0200500 
Hwang and Chen [12] 22.0212 
Marzban and Pirmoradian [23] 22.0230201080 
Proposed method (N = 12) 
ADM 22.02235 
VIM 22.02115 


Example 4. We now consider the following nonlinear time-varying multi- 
delay systems: 


a(t) = x(t — 1)a(t — 2)u(t — 2), 0<t<3, 
x(t) = 1, —2<t<0, (41) 
u(t) = 0, 2 to 0, 
with the cost functional 
3 
J =) (2? (t) + u?(t))dt. (42) 
0 


This optimal control is adopted from [10,34,35]. In order to obtain an accu- 
rate enough suboptimal control law, we applied the proposed algorithm with 
the tolerance bounds ¢ = 1.2 x 107°. Table 6, a comparison is made between 
the value of J obtained by the present method with N = 138, together with 
the value of J reported in the literature. The approximate solution of «(t) 
and u(t), obtained by the proposed method with N = 13 and the results 
of PMP method generated by Gollmann et al. [10] are plotted in Figure 3. 
Therefore, in view of the results, the present method is quite effective. 
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Table 6: The cost functional values for Example 4 


Method Cost functional values 
Wachter et al. [35](600 grid points) 2.763044 
Vanderbei [34](600 grid points) 2.763044 
Gollmann et al. [10](600 grid points) 2.761594156 
Proposed method (N = 13) 
ADM 2.761591012 
VIM 2.761592238 


Figure 3: The suboptimal control and state, when N = 13 for Example 4. 


6 Conclusion 


In this work, the ADM has been successfully applied to find the solution of 
suboptimal control for linear time-varying systems with multiple state and 
control delays, and quadratic cost functional is presented. By using the ADM 
and VIM with the finite-step iteration of algorithm, we obtained a suboptimal 
control law. Some numerical examples have been provided to demonstrate 
the validity and applicability of the proposed method. The method is general 
and yields very accurate results. 
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